First and second order clustering transitions 
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We consider a Hamiltonian system made of TV classical particles moving in two dimensions, cou- 
pled via an infinite-range interaction gauged by a parameter A. This system shows a low energy 
phase with most of the particles trapped in a unique cluster. At higher energy it exhibits a transition 
towards a homogenous phase. For sufficiently strong coupling A an intermediate phase characterized 
by two clusters appears. Depending on the value of A the observed transitions can be either second 
or first order in the canonical ensemble. In the latter case microcanonical results differ dramatically 
from canonical ones. However, a canonical analysis, extended to metastable and unstable states, 
is able to describe the microcanonical equilibrium phase. In particular, a microcanonical negative 
specific heat regime is observed in the proximity of the transition whenever it is canonically discon- 
tinuous. In this regime, microcanonically stable states are shown to correspond to saddles of the 
Helmholtz free energy, located inside the spinodal region. 
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It is well known that thermodynamic quantities de- 
rived within different statistical ensembles should coin- 
cide in the thermodynamic limit. However, this state- 
ment is valid only if the interaction among the parti- 
cles satisfies two conditions : (i) the pairwise interaction 
potential is integrable ( i.e it decays faster than l/r d , 
where d is the space dimension); (ii) the potential en- 
ergy per particle is bounded from below jL|. Whenever 
one of these conditions is violated, ensemble inequiva- 
lence and thermodynamical instabilities can occur. An 
extreme situation is represented by the gravitational po- 
tential for which neither condition is satisfied. Indeed, for 
gravitating systems the usual laws of equilibrium ther- 
modynamics are expected not to hold : one of the most 
striking anomalies is related to the negative values taken 
by the specific heat [Q-Q. Discrepancies between results 
obtained in the microcanonical and canonical ensembles, 
with an associated negative specific heat regime, have 
been observed in the thermodynamic limit for several sys- 
tems with attractive potentials violating either condition 
(i) ||-|7| or condition (ii) only Similar anomalies 

are also present for systems with a finite number of par- 
ticles, e.g. for nuclear multifragmentation |jl0|| , as well 
as for atomic clusters [jll). In all these cases, prelimi- 
nary results suggest that ensemble inequivalence can be 
observed in proximity of a canonically first order phase 
transition |Jl|||. 

In this article, we aim at better clarifying the origin 
of such inequivalence by considering a generalization of 
a previously studied iV-body classical Hamiltonian sys- 



tem with infinite-range attractive interactions |5|]. The 
novelty of the present model consists in the presence 
of a tunable coupling A that allows to change the na- 
ture and the order of the transitions. These are inves- 
tigated analytically within the canonical ensemble and 
numerically via microcanonical molecular dynamics sim- 
ulations. This system shows three different phases : a 
clustered phase [CP\) occuring at low internal energy 
U (temperature T), with most of the particles trapped 
in a unique cluster, a second clustered phase (CP2), ex- 
hibiting two clusters, at intermediate energy (tempera- 
ture) and sufficiently strong coupling A, and a homoge- 
nous phase (HP), with particles uniformly distributed, 
at high values of U (T). The canonical equilibrium solu- 
tions are computed from the lowest lying extrema of the 
Helmholtz free energy and reveal that the system under- 
goes either first or second order transitions depending 
on the value of the coupling constant A. In particular, 
we have focused our attention on first order transitions 
separating the ordered phase CP\ and the homogeneous 
one HP and between the two clustered phases. In both 
cases canonical and microcanonical equilibrium predic- 
tions differ dramatically near the transition, revealing a 
negative specific heat regime within the microcanonical 
ensemble. No discrepancy between microcanonical and 
canonical results appear at continuous phase transitions, 
at least for what concerns the temperature-energy equi- 
librium relation @. 

Our results demonstrate that, irrespectively of the na- 
ture of the phases involved in the transition, the micro- 
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canonical negative specific heat regime can be well re- 
produced even within the canonical ensemble. Once not 
only the absolute minima of the Helmholtz free energy 
are taken in account, but also the relative extrema corre- 
sponding to canonically metastable and unstable states. 

The model we consider is a classical N-body Hamilto- 
nian system defined on a two-dimensional periodic cell. 
The interparticle potential is infinite ranged and all the 
particles are identical with unitary mass. The Hamil- 
tonian of the model is Ha = K + Va where K = 
is the kinetic energy and the potential 



2^1=1 y 2 

energy reads as 



—Acoa(xi - x J )coa(y l - yj)] , 



(1) 



representing the co- 
.,,i,Py,i) the conjugated 



with (xj,2/j) S] — tt : tt] x] — tt : tt 
ordinates of the i-particle and (p. 
momenta. For A = 0, the two spatial directions x and y 
are uncoupled and the Hamiltonian reduces to the sum of 
two independent one-dimensional mean-field models [ p~5[ . 
In this shown in , a second order phase tran- 

sition appears, both in the microcanonical and in the 
canonical ensemble, connecting a single clustered phase 
(CPi), at sufficiently low specific energy U = H/N, to a 
homogeneous phase (HP) at high energy. For non zero 
values of A, the two spatial directions are coupled. Pre- 
vious investigations were limited to the value A = 1 |j| 
and a transition was also observed from a GP\ to a HP 
phase. This transition is first-order in the canonical en- 
semble, while microcanonical simulations are compatible 
with a continuous transition associated with a negative 
specific heat regime. Both for A = and A = 1, at low 
energies, all particles arc trapped in a cluster, while, for 
sufficiently high energies, they are uniformly distributed 
in the cell. 

In order to better investigate the origin of ensemble 
inequivalence within a unique framework, we have intro- 
duced model (|l|) which allows, by continuously varying 
parameter A, to pass from a situation where the micro- 
canonical and canonical results coincide (A = 0) to a 
situation where the two ensembles disagree over a finite 
energy range near the transition. As we will show, this 
model is indeed richer than expected because it reveals 
also more complicated transitions than those previously 
studied in @@. 

Due to the long-range nature of the interactions, 
the collective behaviour of the particles can be de- 
scribed in terms of the following mean-field vectors: 
M z = ((cos(z)) N , (am(z)} N ) = M z (cos(^), sin(<^)) 
where <p z G [0, 7r/2] and z — x or y; P ' z = 
((cos(z)) JV , (sin(z)) w ) = P z (cos(V> z ),sin(-0 z )) where 
ijj z £ [0, 7r/2] and z = x ± y. The average over all the 
particles is indicated by (■■■) N . It can be shown that on 



average M x w My w M and P x + y ~ Px-y ~ P (for 
more details see ||). Therefore, the potential energy 
can be rewritten, in the mean- field limit N — » oo, as 
V A = [2 + A - 2M 2 - AP 2 }/2. 

For U « (or equivalently at low temperature), the 
system described by model (|l|) is in the CP\, particles 
have all the same location in a single point-like clus- 
ter and M w P rs 1, whereas at large enough energy 
(temperature) the system is in the HP and M sw P = 
0(1/Vn). For sufficiently high values of A > A 2 — 3.5, 
a third intermediate phase CP2, exhibiting two clusters, 
appears. In this phase, due to the symmetric location of 
the two clusters M - 0(1/ VF) and P - 0(1). 




FIG. 1. Phase diagrams of model (|lj) reporting the tran- 
sition temperatures versus the coupling parameter A (a) and 
the corresponding specific energy U = H/N versus A (b). 
The solid lines indicate the canonical transition lines and the 
dots the points where the nature of the transitions change. 
A\, A2 and A3 are the threshold coupling constants that de- 
termine the transition scenario I — > IV displayed above the 
graphs. The grey shaded area in (b) indicates the domain 
where two phases coexist. Inside this area, the dashed curves 
are the two spinodal lines. 



In the mean-field limit, the equilibrium properties of 
model (Q) can be derived analytically within the canon- 
ical ensemble following the approach of Ref . || . In par- 
ticular, the Helmholtz free energy reads as 



F(M,P;T,A) 



= T(M 2 + P 2 ) -ln[T G{M,P;A)} (2) 



with G = Jo I (M + V2APcos(s))exp(Mco$(s))ds, 
where Iq is the modified Besscl function of zero order. 
Since we are interested also in metastable and unstable 
states we will not restrict ourselves to the study of the 
lowest-lying minima of F, but we will keep track of all 
the other extrema. Due to the P — > —P symmetry of 
F, we can limit our analysis to the P > half-plane. 
The transition lines, obtained by considering the abso- 
lute minima of F, are reported as solid lines in Fig. [I]. 
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The minima can be easily associated to the three ob- 
served phases, since HP will correspond to M = P = 0, 
CPi to \M\ > 0,\P\ > and CP 2 to M = and \P\ > 0. 



FIG. 2. Canonical temperature-energy relation in the co- 
existence region for A=l (a) and A=4 (b). Lines indicate 
canonical analytical results, while circles correspond to mi- 
crocanonical simulations. Solid thick lines are equilibrium re- 
sults, solid thin lines metastable states and dashed thin lines 
unstable states. The straight dashed thick line is the Maxwell 
construction. Figure (a) refers to a first order transition from 
CPi to HP, (b) to discontinuous transition connecting the 
two clustered phases. In (b) the second order transition from 
CP2 to HP is also shown. The microcanonical results have 
been obtained via molecular dynamics simulations of model 
(§) with N = 5, 000 particles (is) . 



Let us describe the observed phase transitions within 
the canonical ensemble with the help of Fig. |l|: the line 
referred to as Tm in the inset (a) indicates the transi- 
tion temperatures where M vanishes and where the phase 
CP\ looses its stability; while the line Tp is where P — > 
and the CP 2 leaves place to the HP. The phase diagram 
U versus A reported in the inset (b) gives clearer hints 
for what concerns first order transitions, indicating the 
corresponding energy jumps (latent heats). Depending 
on the value of A, four different scenarios can be distin- 
guished. (I): [0 < A < A\ = 2/5] - in this case one 
observes a continuous transition from a CPi to a HP, 
the critical line is located at Tm = 1/2 (Um = 3/2 + A). 
Along this line, at A\ = 2/5, there is a tricritical point 
: here the transition becomes discontinuous. (II) : 
[Ai < A < A 2 w 3.5] - The transition between CPi 
and HP is first order in all this range of parameters with 
a finite energy jump (grey shaded area in Fig. 0(b)). 
(Ill) : [A 2 < A < A 3 « 5.7] - In this region the third 
phase begins to play a role and two successive transitions 
are observed: first CP\ disappears at Tm via a first order 
transition that gives rise to the biclustered phase CP 2 , 
which ends up in the HP due to a continuous transi- 
tion. The critical line associated to this last transition is 
T P = A/ A (Up = 3A/4 + 1). (IV) : [A > A 3 ] - Finally, 
region (IV) differs from region (III) just for the nature of 
the transition connecting the two clustered phases, that 
becomes second order. 
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Let us now concentrate on first order phase transi- 
tions, in particular we will compare microcanonical re- 
sults obtained via molecular dynamics simulations Eg ] 
with canonical results. Outside the grey area in Fig. 0(b) 
the microcanical and canonical results coincide every- 
where, but inside such region strong discrepancies are 
observed. This can be clearly seen in Fig. |[ where T is 
plotted as a function of U near the discontinuous tran- 
sition for two different values of A (namely, A = 1 and 
A = 4). Within the canonical ensemble, first order phase 
transitions occur at a given temperature and the coexis- 
tence region is bridged with a horizontal line (T = Tm in 
Fig. |(a) and T = T P in |(b)), according to the Maxwell 
construction. 

However, this region is accessible via microcanonical 
simulations and the corresponding equilibrium results are 
shown as circles in Fig. 0. Within this region these re- 
sults disagree with the canonical ones, the transition is 
now continuous and a negative specific heat region ap- 
pears. It is remarkable that both in the transition con- 
necting the clustered phase CP\ to the HP as well as 
in that connecting the two clustered phases the essential 
features of the transition are the same. 

We will now show that the microcanonical results can 
be obtained even within the canonical ensemble by con- 
sidering also the relative minima and the saddles of F. 
A schematic sketch of the relevant extrema of F in the 
(M, P)-plane close to the two transitions discussed above 
is reported in Fig. ||. Let us first consider the case A = 1 
for 1/2 < T < 0.551, that is reported in Fig. |(a). For 
T < 1/2 the free energy has a unique minimum corre- 
sponding to the CP\ phase, while at T = 1/2 a second 
relative minimum emerges at M = P = associated to a 
metastable HP state. Exactly at the same temperature 
also two symmetric saddles appear in the [M, i- > )-plane 
inbetween the CP\ and the HP minima. For increasing 
temperature the depth of the CP\ minimum decreases 
while that of the HP minimum increases and, simulta- 
neously, the saddles and the CP± minima approach each 
other. At T = T M ^ 0.54 the CPi and HP mimima 
reach exactly the same depth and this singles out the 
canonical transition temperature. At T = 0.551 the sad- 
dles and the CP± minima finally merge and disappear. If 
the energy and temperature denoting these unstable and 
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metastable states are reported in a graph together with 
the equilibrium solutions, one observes that they con- 
nect continuously the CP\ equilibria to the HP ones (see 
Fig. ||). Moreover, as shown again in Fig. ||(a), the micro- 
canonical results essentially coincide with the metastable 
and unstable canonical solutions in the coexistence re- 
gion. In particular, the negative specific heat regime is 
associated to the saddles bridging the metastable min- 
ima. The limits of existence of these unstable states 
are tipically referred to as spinodals and are reported 
in Fig. |l| as dashed lines fnj. A similar behaviour is 
observed for the CP\ to CP2 transition, in this case for 
T < 0.93 the unique minimum of F is again associated to 
the CP\ and at T = 0.93 two symmetric minima appear 
on the M = axis with \P\ 7^ 0, and these metastable 
states are clearly associated to the CP2 ■ Also two saddle 
points appear in the F profile separating the CP\ min- 
ima from the metastable CP2 minimum in the positive 
P semi-plane. For T > Tp ~ 0.95 the CP2 minima be- 
come the stable ones. For increasing temperatures the 
CP\ minima and the saddles approach each other and fi- 
nally vanish at the spinodal (located at T ~ 0.965). For 
higher temperatures the CP^s approach symmetrically 
the origin, where a saddle (not involved in the first order 
transition) is present corresponding to the HP unstable 
phase. At T = Tm — 1 the two minima merge with 
the saddle at the origin, giving rise to a unique mini- 
mum for F corresponding to a stable HP. This second 
transition is clearly continuous. Also in this case, the 
microcanonical continuous transition is well reproduced 
by the metastable and unstable states involved in the 
discontinuous transition (see Fig. 0(b)). 
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FIG. 3. Schematic representation of the extrema of the 
free energy for a fixed temperature in the intervals reported 
in the figures for A = 1 (a) and A = 4 (b). Circles denote the 
positions of the minima, while the crosses represent the sad- 
dles. The motion of the extrema for increasing temperatures 
is indicated by the arrows. 

In conclusion, we have shown that a study of the ab- 
solute and relative extrema of the Helmholtz free energy 
is sufficient to provide a complete microcanonical and 
canonical description of the equilibrium behaviour of a 
iV-body Hamiltonian with infinite-range attractive inter- 
action. Depending on the value of the tuning parameter 
A canonically first or second order transitions are ob- 



served. The comparison of the canonical results with the 
microcanonical ones shows that ensemble inequivalence 
occurs in proximity of canonically discontinuous transi- 
tions. Irrespective of the phases involved in the transi- 
tions a negative specific heat regime has been observed. 
In particular, this regime is always associated to the ex- 
istence of a spinodal region for the canonical solutions. 
For energy values within this region microcanonically sta- 
ble states corresponds to saddles of the free energy. We 
have validated this scenario by considering two differ- 
ent discontinuous transitions connecting the one cluster 
phase either to the homogeneous one or to the two cluster 
phase. 

We believe that the results reported in this article are 
not limited to systems with infinite-range interactions, 
but they should be applicable also to systems with power- 
law decaying potentials (as shown in |18[) and to finite 
systems with short-ranged forces (see [12|). 
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